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WEIGHT  STUDY  OF  PARTIALLY  SEGMENTED  DIRECT-CONDENSING  RADIATORS 

FOR  LARGE  SPACE  POWER  SYSTEMS 


by  Roy  L.  Johnson 
Lewis  Research  Center 

SUMMARY 

An  analytical  investigation  has  been  performed  to  provide  an  insight  into  the  weight 
increase  of  a  direct -condensing  radiator  for  a  potassium  Rankine  cycle  space  power 
system  as  the  system  size  was  increased  from  1  to  10  megawatts  electric  while  holding 
constant  the  component  efficiencies  and  all  cycle  temperatures  and  pressures. 

The  largest  single  contributor  to  the  radiator  weight  is  the  armor  necessary  to  pro¬ 
tect  fluid  passages  from  meteoroid  puncture.  The  two  most  important  factors  affecting 
armor  thickness  are  the  meteoroid  population  estimate  and  the  radiator  design  survival 
probability.  Design  survival  probabilities  of  0.  9,  0.  95,  and  0.  99  were  used  with  a  high 
meteoroid  population  estimate,  Whipple’s  (1961),  and  a  low  estimate,  Watson's  (1956), 
for  1,  10,  20,  and  100  radiator  segments. 

The  radiator  weight  is  sharply  influenced  by  the  meteoroid  flux  estimate  and  the  de¬ 
sign  survival  probability.  For  the  unsegmented  case,  using  Watson's  (1956)  flux  esti¬ 
mate,  the  radiator  weight  ranged  from  1.  5  to  3  pounds  per  kilowatt  electric  for  the  1-  to 
10-megawatt  power  range  for  a  0.  9  radiator  design  survival  probability,  and  2.  5  to 
6  pounds  per  kilowatt  for  a  0.  99  design  survival  probability.  When  Whipple's  (1961)  flux 
estimate  was  used,  the  radiator  weight  ranged  from  7  to  17  pounds  per  kilowatt  for  a 
0.  9  design  survival  probability  and  from  16  to  42  pounds  per  kilowatt  for  a  0.  99  design 
survival  probability . 

Segmentation  of  the  radiator  provided  significant  weight  savings  over  the  unseg¬ 
mented  radiator,  especially  when  high  survival  probabilities  were  demanded.  Generally 
the  weight  advantages  of  segmentation  increased  with  increases  in  system  power  level, 
design  survival  probability,  number  of  segments,  and  severity  of  the  meteoroid  flux  es¬ 
timate.  When  Whipple's  (1961)  meteoroid  flux  estimate  and  a  0.  99  design  survival  prob¬ 
ability  for  a  10-megawatt  system  were  used,  the  radiator  with  100  segments  was  about 
one-fourth  the  weight  of  an  unsegmented  radiator. 

With  segmentation  the  ratio  N  /N  of  the  design  value  of  surviving  number  of  seg- 
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ments  to  the  total  number  of  segments  that  yielded  minimum  weight  was  found  to  be  0,8 
or  higher.  In  this  study,  however,  the  manifolds  were  not  segmented.  Had  the  manifold 
been  segmented  N  /N  would  tend  to  be  lower  and  conceivably  might  approach  the  0.  75 
previously  reported  in  the  literature. 


INTRODUCTION 


)  The  heat  rejection  aspects  of  electric  generating  systems  for  space  application  have 
received  considerable  attention  in  recent  years.  It  is  recognized  that  a  large  fraction  of 
the  power  system  weight  would  be  taken  up  by  the  waste  heat  radiator  if  the  conventional 
fin  and  tube  type  were  used.  Several  radiator  studies  have  been  reported  in  the  literature 
for  systems  unto  several  megawatts  of  generated  powerj  A  partial  list  includes  refer¬ 
ences  1  to  6.  [it  would  be  difficult  to  relate  the  results  of  these  studies  to  one  another 
since  each  study  used  different  temperature  levels,  materials,  and  criteria  for  calcula¬ 
ting  protective  armor  thicknesses. 

The  analytical  study  presented  herein  is  an  attempt  to  determine  weight  trends  that 
might  be  expected  for  both  segmented  and  unsegmented  direct-condensing  fin  and  tube 
radiators  for  a  Rankine  cycle  turboelectric  system  for  the  range  of  power  levels  from  1 
to  10  megawatts  electric.  No  allowance  was  made^or  the  weight  of  any  necessary  radia¬ 
tor  support  structure.  The  working  fluid  was  potassium  and  the  materials  used  were 
beryllium  for  the  fin^  and  tube  ari^or  and  columbium  for  the  tube  li^r .  When  high  and 
low  metec^id  flux  estimates  were  used,  weights  obtained  when  segmentation  of  the  tube 
panels  was  employed  were  compared  to  the  unsegmented  radiator  weights.  Of  interest 
also,  was  the  optimum  ratio  of  design  value  of  surviving  segments  to  the  total  number  of 
segments^  Reference  7  had  determined  the  optimum  ratio  of  surviving  to  total  number  of 
segments  for  an  armored  radiator  without  fins  and  manifolds. 

pfhe  study  objective  was  to  glean  trends  from  many  detailed  radiator  designs  that 
could  be  used  to  obtain  rough  estimates  of  the  effect  of  power  level  and  segmentation  on 
radiator  weight  without  having  to  make  detailed  calculations.  Therefore,  for  compara- 
,  identical  cycle  conditions  were  maintained  throughout  the  power  range 
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tive  purpose 
considered.! 


PROCEDURES  AND  ASSUMPTIONS 

The  procedure  and  the  equations  used  are  covered  in  some  detail  in  appendix  B.  The 
development  of  the  equations  used  to  calculate  such  items  as  fin  thickness,  fin  length, 
manifold  weight  and  surface  area,  tube  and  manifold  pressure  drops,  and  fin  and  tube 
angle  factors  are  shown  in  appendixes  C  to  H. 
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It  was  assumed  that  saturated  vapor  at  1960°  F  entered  the  turbine,  the  turbine  ef¬ 
ficiency  was  0.8,  the  generator  efficiency  was  0.9,  and  a  turbine  exit  and  condensing 
temperature  of  1400°  F  had  the  corresponding  vapor  pressure  of  about  15  psia.  The  tur¬ 
bine  temperature  ratio  was  taken  from  reference  8  for  minimum  radiator  area  for  the 
1400°  F  condensing  temperature.  Selection  of  cycle  temperatures  and  efficiencies  of  the 
turbine  and  generator  dictates  the  turbine  exit  quality,  fluid  mass  flow,  and  radiator  heat 
load.  The  calculated  turbine  exit  quality  was  0. 85,  the  fluid  mass  flow  was  6. 88  pounds 
per  second  per  megawatt  electric,  and  the  radiator  heat  load,  including  100°  F  subcool¬ 
ing,  was  l.OOlTxlo"^  Btu’s  per  hour  per  megawatt  electric. 

The  following  is  a  brief  discussion  of  other  assumptions  and  restrictions  that  in¬ 
fluence  the  study  results. 


Materials 

The  only  materials  considered  were  beryllium  for  the  fins  and  meteoroid  protective 
armor,  and  a  0.  03-inch  columbium  liner  for  the  fluid  passages. 


Radiator  Panel  Configuration 

Two  configurations  were  used,  both  radiating  from  both  sides.  One  was  a  two-panel 
configuration  (fig.  1(a))  with  a  vapor  manifold  at  one  end  and  a  liquid  manifold  at  the 
other  end.  The  other  configuration  (fig.  1(b))  with  four  panels  had  a  central  vapor  mani¬ 
fold  and  liquid  manifolds  at  the  outer  extremities.  All  manifolds  were  tapered. 


(a)  Two-pane!  configuration. 


(b)  Four-panel  configuration. 


Figure  1.  -  Radiator  panel  arrangement. 
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Angle  Factors 


In  the  radiant  heat-transfer  treatment,  interradiation  and  interreflection  between 
surfaces  were  not  considered.  The  angle  factors  for  both  the  fin  and  tube  took  into  ac¬ 
count  only  the  direct  radiation  to  space.  The  error  introduced  by  such  an  assumption  has 
been  discussed  in  reference  9,  and  the  maximum  error  for  the  results  reported  herein 
was  approximately  8  percent.  That  is,  the  calculated  radiant  heat  transferred  from  the 
fins  and  tubes  was  less  than  the  actual  amount  by  a  maximum  of  8  percent. 


Fins  and  Tubes 

Only  constant  thickness  fins  and  constant  diameter  tubes  were  considered  even 
though  tapered  fins  and  tubes  may  have  offered  some  savings  in  weight.  Reference  10 
shows  a  maximum  weight  savings  of  10  percent  with  tapered  fins,  but  did  not  include 
manifold  weights.  Had  the  manifold  weights  been  included,  the  weight  savings  would 
have  been  even  less. 


Factors  Affecting  Armor  Thickness 

Radiator  weight  is  closely  tied  to  the  amount  of  armor  necessary  to  protect  against 
meteoroid  damage.  The  uncertainty  of  the  numbers  of  particles  that  would  be  encountered 
and  of  the  damage  that  would  be  inflicted  by  any  one  meteoroid  particle  makes  a  calcula¬ 
tion  of  armor  thickness  somewhat  speculative.  Several  estimates  of  the  meteoroid  popu¬ 
lation  are  shown  graphically  in  reference  11  and  differ  by  orders  of  magnitude.  The  es¬ 
timates  for  the  density  of  meteoroids  vary  from  0. 05  to  7.9  grams  per  cubic  centimeter 
(ref.  12).  Due  to  these  wide  differences  in  meteoroid  flux  and  mass  estimates,  two  sets 
of  meteoroid  data  from  reference  11  were  used  in  order  to  encompass  a  range.  Arbi¬ 
trarily  selected  were  Whipple's  (1961)  and  Watson's  (1956)  meteoroid  flux  estimates.  The 
meteoroid  flux  and  mass  estimates  and  the  other  important  inputs  used  in  the  calculation 
of  armor  thickness  are  the  following: 

(1)  Meteoroid  flux. 

(a)  High  estimate,  Whipple's  (1961);  =  9.27x10”"^  particles/ (ft^)(yr) 

(eq.  for  curve  in  ref.  11).  The  meteoroid  density  recommended  for  use  with  this 
flux  was  Pp  =  2.  7  grams  per  cubic  centimeter. 

(b)  Low  estimate,  Watson's  (1956):  =  1. 19x10”®  mp^  particles/ (ft^)(yr)  (eq. 

for  curve  in  ref.  11).  A  meteoroid  density  of  Pp  =  0.  5  gram  per  cubic  centimeter 
was  arbitrarily  assumed.  The  armor  thickness  is  only  weakly  influenced  by  particle 
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density.  The  meteoroid  flux  used  is  the  dominant  factor  in  influencing  the  calculated 
armor  thickness  deemed  necessary. 

(2)  Penetration  criterion.  The  model  proposed  in  reference  11  was  used.  The  equa¬ 
tion  is  shown  as  equation  (B7 )  of  appendix  B . 

(3)  Survival  probability,  0.9,  0.95,  and  0.99. 

(4)  Mission  duration,  2  years . 


Manifold  Fluid  Velocity 

Tapered  manifolds  maintained  constant  fluid  velocities.  The  fluid  velocity  was 
chosen,  then  the  manifold  internal  diameter  was  calculated  to  maintain  the  velocity 
picked.  For  most  of  the  study  the  fluid  velocities  used  were  200  feet  per  second  in  the 
vapor  manifold  and  2  feet  per  second  in  the  liquid.  If  the  sum  of  the  calculated  manifold 
pressure  drops  exceeded  0.  25  psi,  the  velocities  were  reduced  to  150  feet  per  second 
for  the  vapor  and  1  foot  per  second  for  the  liquid.  If  the  sum  of  the  pressure  drops  still 
exceeded  0. 25  psi,  the  vapor  velocity  was  further  reduced  to  100  feet  per  second. 


Other  Assumptions 

Radiant  input.  -  The  effective  sink  temperature  was  assumed  to  be  0°  R. 

Working  fluid  temperature.  -  The  condensing  working  fluid  was  assumed  to  be  iso¬ 
thermal  everywhere.  The  fluid  stagnation  temperature  was  used  to  obtain  the  fluid  ther¬ 
modynamic  properties. 

Internal  tube  wall  temperature.  -  The  tube  internal  wall  temperature  was  assumed 
to  be  equal  to  the  fluid  stagnation  temperature. 

Subcooling.  -  The  condensed  liquid  was  subcooled  100°  F. 

Emissivity.  -  The  external  radiating  surface  emissivity  was  e  =  0.  9. 


PARAMETERS  FIXED  AS  A  RESULT  OF  PRELIMINARY  CALCULATIONS 

Preliminary  calculations  were  made  in  an  attempt  to  reduce  the  number  of  variables 
to  be  used  throughout  the  study.  A  discussion  of  these  preliminary  findings  and  reasons 
for  picking  the  values  selected  follows. 


Generalized  Length  Parameter 

In  reference  13  the  value  for  the  generalized  fin  length  parameter  that  yielded 
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.6  .8  1.0  1.2 
Generalized  length  parameter,  Lg 

Figure  2.  -  Effect  of  generalized  length 
parameter  on  radiator  weight  and  area. 
Whipple's  flux  (1961);  system  power  out¬ 
put.  10  megawatts;  probability  of  re¬ 
ceiving  no  meteoroid  puncture  during 
mission  duration,  0.9;  individual  tube 
length,  40  feet;  tube  static  pressure 
drop,  1  psi;  two-panel  configuration. 


Figure  3.  -  Specific  weight  and  ratio  of  fin  span  to  tube  outside  diameter  as 
function  of  number  of  tubes.  Watson's  flux  (1956);  system  power  output, 

1  megawatt  electric;  probability  of  receiving  no  meteoroid  puncture  during 
mission  duration,  0.  9;  individual  tube  length,  15  feet;  tube  static  pres¬ 
sure  drop,  1  psi;  generalized  length  parameter,  1;  two-panel  configuration. 


the  minimum  volume  fin  was  approximately  0.  92  for  the  case  where  the  sink  temperature 
was  0°  R.  The  value  of  that  would  yield  the  lightest  weight  radiator  when  tubes  and 
manifolds  were  included  was  not  obvious.  Values  for  L  of  0.7,  0.92,  1.0,  and  1. 1 

o 

were  used  to  make  some  preliminary  radiator  weight  and  area  calculations.  For  each 
value  of  Lg  the  number  of  tubes  was  varied  until  a  minimized  weight  and  corresponding 
tube  diameter  were  obtained,  while  everything  else  such  as  individual  tube  length,  radi¬ 
ator  heat  load,  and  tube  pressure  drop  were  held  constant.  Figure  2  is  typical  of  the  re¬ 
sult  of  plotting  the  percent  deviation  from  minimum  weight  against  Lg.  For  these  pre¬ 
liminary  calculations,  the  minimum  weight  usually  occurred  at  a  value  of  Lg  close  to 
unity.  Figures  11  and  12  of  reference  10  showed  combined  fin  and  tube  heat  rejection 
per  unit  weight,  plotted  against  the  ratio  of  half  fin  length  to  tube  outer  radius,  to  be 
maximum  for  a  conductance  parameter  of  about  1.  The  conductance  parameter  of  refer¬ 
ence  10  is  identical  to  Lg^e.  Using  an  emissivity  e  of  0.  9  the  conductance  parameter 
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of  unity  (of  ref.  10)  corresponds  to  a  value  for  Lg  of  about  0.  95.  In  reference  10,  the 
particular  value  of  which  the  maximum  heat  rejection  occurred  was  differ¬ 

ent  for  each  case  shown. 

While  the  approach  taken  in  the  study  reported  herein  was  somewhat  different  from 
that  of  reference  10,  changing  the  number  of  fixed  length  tubes,  while  maintaining  a  con¬ 
stant  L  ,  changes  the  percentages  of  the  total  heat  load  rejected  by  the  tubes  and  by  the 
fins.  This  in  essence  is  the  same  as  changing  the  ratio  This  can  be  seen 

from  figure  3  where  specific  weight  and  21  are  plotted,  for  one  example,  against 

the  number  of  fixed  length  tubes.  Decreasing  increases  D.  and  but 

also  increases.  Thus  for  finding  minimum  radiator  weight  (constant  Lg),  changing  N^. 
is  equivalent  to  changing  All  of  the  subsequent  radiator  weights  were  obtained 

by  using  a  value  of  L  =1.  While  a  value  close  to  unity  provided  the  minimum  weight, 

it  should  not  be  overlooked  that  by  using  a  lower  value  of  Lg  the  projected  radiator  area 
may  be  decreased  somewhat  for  only  a  slight  radiator  weight  penalty.  For  example,  re¬ 
turning  to  figure  2,  when  a  value  of  Lg  =  0. 7  is  used,  the  area  could  be  reduced  by  per¬ 
haps  9  percent  with  a  weight  increase  of  about  2  percent.  For  some  applications  it  might 
be  desirable  to  accept  the  higher  weight  to  get  the  smaller  area. 


Pressure  Drop 

The  pressure  drop  across  the  radiator  affects  its  weight  in  two  ways.  The  size  of 
the  fluid  passages  is  a  function  of  the  pressure  drop,  a  high  pressure  drop  permitting 
small  passages.  High  pressure  drop,  however,  reduces  the  average  condensing  temper¬ 
ature  as  well  as  adversely  affecting  the  weight  of  the  whole  power  system  by  depressing 
the  cycle  efficiency;  this  effect  has  been  ignored  in  this  analysis.  Furthermore,  a  high 
radiator  pressure  drop  results  in  a  low  condensate  pump  inlet  pressure  which  may  make 
necessary  additional  subcooling  to  prevent  pump  cavitation  problems.  Figure  4  shows  a 
weight  proportionality  factor,  derived  by  dividing  radiator  weight  by  the  weight  of  an 
arbitrarily  chosen  radiator,  for  three  tube  pressure  drops  plotted  against  the  number  of 
tubes  for  one  power  level  as  an  illustrative  example.  The  minimum  weights  for  each 
tube  length  were  about  the  same  for  the  1-  and  2-psi  tube  static  pressure  drops  and 
higher  for  the  5-psi  pressure  drop.  From  figure  4  a  tube  static  pressure  drop  of  1  psi 
was  chosen  for  the  remainder  of  the  study.  This  may  not  be  the  optimum  pressure  drop, 
but  the  trends  being  sought  should  not  be  influenced  significantly .  Another  0. 25  psi  max¬ 
imum  was  allowed  for  the  sum  of  the  pressure  drops  in  the  vapor  and  liquid  manifold. 
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Figure  4.  -  Effect  of  pressure  drop,  tube  length,  and  number  of  tubes  on 
weight.  Meteoroid  flux,  system  power  output,  probability  of  receiving 
meteoroid  puncture  during  mission  duration,  generalized  length  param¬ 
eter,  and  number  of  panels  are  constant. 


DISCUSSION  OF  RESULTS 

Radiator  calculations  were  made  at  four  power  levels:  1,  2,  5,  and  10  megawatts 
electric.  For  each  power  level  several  tube  lengths  were  considered,  and  for  each  tube 
length  specific  radiator  weight  was  plotted  against  the  number  of  tubes.  This  was  done 
for  the  unsegmented  radiator  for  all  combinations  of  the  three  design  survival  probabil¬ 
ities,  two  meteoroid  flux  estimates,  and  the  two  panel  configurations  considered.  These 
plots  are  too  numerous  to  include,  but  figure  5  is  included  as  a  representative  example. 
Tube  length  and  the  number  of  tubes  were  varied  until  the  combination  that  yielded  mini¬ 
mum  weight  was  obtained  for  each  power  level.  In  some  instances  the  minimum  weight 
for  a  particular  tube  length  could  not  be  obtained  because  the  calculations  were  stopped 
when  the  fluid  velocity  at  the  tube  inlet  reached  1300  feet  per  second.  Reference  14  in¬ 
dicates  that  the  vapor  velocity  for  critical  two-phase  flow  (steam  water  mixtures)  is 
lower  than  the  sonic  velocity  for  the  vapor.  The  equilibrium  sonic  velocity  for  potassium 
saturated  vapor  (ref.  15)  for  the  1400°  F  condensing  temperature  is  about  1450  feet  per 
second.  The  1300  feet  per  second  limit  was  chosen  arbitrarily  to  stay  below  the  vapor 
equilibrium  sonic  velocity  by  a  safe  margin. 
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Figure  5.  -  Effect  of  tube  length  and  number  of  tubes  on  radiator  specific  weight. 
Whipple's  flux  (1961);  unsegmented;  probability  of  receiving  no  meteoroid  punc¬ 
ture  during  mission  duration,  0.95;  tube  static  pressure  drop,  1  psi;  gener¬ 
alized  length  parameter,  1.0;  two-panel  configuration. 


Effect  of  Number  of  Tubes  and  Tube  Length  on  Tube  Internal  Diameter 


The  internal  diameter  of  the  tube  for  any  particular  radiator  may  be  of  interest  be¬ 
cause  of  its  influence  on  the  area  vulnerable  to  meteoroid  puncture  as  well  as  heat- 
transfer  considerations.  As  previously  mentioned,  a  constant  tube  static  pressure  drop 
of  1±0. 05  psi  was  used  for  the  entire  study.  The  internal  tube  diameter  and  approximate 
fluid  velocity  at  tube  inlet  for  any  radiator  shown  in  figure  5  can  be  obtained  from  fig¬ 
ure  6  (p.  10).  It  is  seen  from  figure  5  that  for  any  given  power  level  increasing  the  tube 
length  decreases  the  number  of  tubes  at  which  the  minimum  weight  occurs  and  that  indi¬ 
vidual  internal  tube  diameter  increases  (see  fig.  6).  Figure  6  also  shows  that  for  a 
given  tube  length  the  internal  tube  diameter  is  relatively  insensitive  to  changes  in  the 
number  of  tubes,  even  though  the  fluid  velocity  at  tube  inlet  varies  considerably .  This 
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Figure  6.  -  Internal  tube  diameter  as  function  of  tube  length  and  number  of  tubes.  Tube  static  pressure  drop,  1  psT 


indicates  that  as  inlet  velocity  increases,  increases  in  friction  pressure  drop  are  offset 
by  increases  in  the  momentum  pressure  rise. 


Planform  Geometry 

The  ratio  of  the  width  of  the  radiator  to  the  length  in  the  vehicle  axial  direction  may 
be  important  from  a  vehicle  integration  standpoint.  Flexibility  in  choosing  either  dimen¬ 
sion  without  having  to  pay  a  severe  weight  penalty  would  be  desirable.  Plots  of  specific 
weight  against  the  ratio  of  total  width  to  radiator  length  (fig.  7)  indicate  that  some  lati¬ 
tude  in  the  choice  of  planform  layout  does  exist.  The  minimum  weight  points  from  plots 
similar  to  figure  5  were  plotted  to  show  specific  weight  as  a  function  of  the  ratio  of  radi¬ 
ator  width  to  length  in  figure  7.  Each  curve  represents  the  locus  of  the  minimum  weight 
points. 

For  some  of  the  curves  the  specific  weight  was  not  very  sensitive  to  the  width  to 
length  ratio.  To  cite  one  example,  in  figure  7(a)  for  the  four -panel  configuration  for 
2  megawatts  electric  and  0.  9  design  survival  probability  using  Whipple’s  (1961)  meteo¬ 
roid  flux  estimate,  the  width  to  length  ratio  ranged  from  0.  4  to  10  with  a  specific  weight 
variation  of  only  about  10  percent.  Generally  it  can  be  said  that  as  the  survival  probabil¬ 
ity  is  increased  or  the  meteoroid  flux  estimate  is  increased,  increasing  the  meteoroid 
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protection  problem,  sensitivity  to  changes  in  width  to  length  ratio  decreases  and  the 
minimum  weight  panel  tends  to  approach  a  square. 


Influence  of  Power  on  Weight 

The  minimum  weights  for  each  combination  of  power  level,  survival  probability, 
meteoroid  flux,  and  number  of  panels  were  taken  from  plots  similar  to  figure  5  (p.  9) 
and  used  as  the  basis  of  the  discussion  in  the  section  that  follows. 

When  the  minimum  radiator  weights  at  each  of  the  four  system  power  levels  were 
plotted  against  the  system  electrical  power,  the  points  on  a  log-log  plot  fell  in  a  straight 
line,  which  suggested  an  expression  of  the  form 

W  oc  (1) 


System  power  output,  Pg,  MWe 

(a)  Two-panel  configuration.  (b)  Four-panel  configuration. 


Figure  8.  -  Effect  of  meteoroid  flux  estimate,  survival  probability,  and  system  power  level  on  radiator  specific  weight. 


as  used  in  reference  5.  (All  symbols  are  defined  in  appendix  A. )  If  the  value  of  the  con¬ 
stant  exponent  b  is  known,  radiator  weights  can  be  scaled  from  one  power  level  to 
another  with  some  degree  of  confidence. 

Equation  (1)  was  rearranged  to 


which  is  the  form  of  the  curves  plotted  in  figure  8  for  the  unsegmented  radiator.  The  use 
of  Whipple's  (1961)  meteoroid  flux  estimate,  for  the  0.  9  radiator  survival  probability, 
yields  radiator  specific  weights  of  7  pounds  per  kilowatt  electric  for  a  1-megawatt  system 
up  to  18  pounds  per  kilowatt  electric  for  a  10-megawatt  system.  For  the  0.  9  survival 
probability,  using  Watson’s  (1956)  meteoroid  flux  estimate,  the  radiator  specific  weight 
ranges  from  about  1.3  pounds  mass  per  kilowatt  electric  for  a  1-megawatt  system  to 
about  3  pounds  mass  per  kilowatt  electric  for  a  10-megawatt  system.  For  the  higher 
0.  99  survival  probability,  the  specific  weights  were  heavier  than  the  specific  weights  as¬ 
sociated  with  the  0.  9  survival  probability  by  a  factor  of  two  or  more. 

The  value  of  b  in  equation  (2)  was  approximately  1.4  when  Whipple's  (1961)  meteo¬ 
roid  flux  estimate  was  used,  and  only  slightly  less  when  Watson's  (1956)  meteoroid  flux 
estimate  was  used.  Increasing  the  design  survival  probability  from  0. 9  to  0.  99  increased 
the  value  of  the  exponent  b  only  about  2  percent.  Also,  the  value  of  b  was  essentially 
the  same  for  both  the  two-  and  four -panel  configurations.  The  value  of  the  exponent  re¬ 
ported  in  reference  5  was  1.  for  a  beryllium  radiator  for  a  heat  load  between  1  and 
10  megawatts  thermal.  Results  plotted  in  reference  6  would  lead  to  a  value  of  less  than 
1. 2  for  a  direct  condensing  radiator  and  slightly  higher  for  a  liquid  filled  radiator,  both 
using  copper  as  the  fin  material. 


Segmenting 

Reference  7  indicates  that  segmenting  the  radiator  to  reduce  the  amount  of  necessary 
protective  armor  should  offer  significant  weight  savings.  Affecting  the  amount  of  pro¬ 
tective  armor  necessary  for  a  given  overall  survival  probability  are  the  number  of  seg¬ 
ments  N  and  the  design  value  for  the  number  of  segments  remaining  unpunctured  at 
the  end  of  the  mission.  Probabilities  S  of  0.  9,  0.  95,  and  0. 99  that  at  least  out  of 
N  segments  would  not  be  punctured  for  the  2-year  mission  duration  were  assumed  to 
calculate  the  armor  thickness  for  varying  while  N  was  held  constant  at  10,  20,  and 
finally  100  segments.  It  should  be  pointed  out  that  the  manifolds  were  not  segmented, 
only  the  tube  banks,  and  that  calculations  were  made  for  the  two-panel  configuration  only. 
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The  equation  for  calculating  the  armor  thickness  is  shown  in  appendix  B.  Furthermore, 
it  was  stipulated  that  all  of  the  radiator  heat  load  was  to  be  rejected  from  segments 
and  the  manifolds. 

Reference  7  showed  a  figure  of  merit,  which  is  proportional  to  weight,  plotted  as  a 

function  of  the  ratio  of  N_  to  N  from  which  it  was  concluded  that  to  obtain  minimum 

s 

radiator  weight  the  optimum  ratio  of  N  to  N  was  approximately  3/4.  Reference  7  did 

o 

not  include  manifolds  and  fins.  An  attempt  was  made  to  compare  some  of  the  results  of 
reference  7  with  the  results  of  this  study  which  does  include  manifolds  and  fins. 

In  reference  7  it  was  assumed  that  the  required  armor  thickness  t^  was  proportional 
to  the  cube  root  of  meteoroid  particle  mass  m  ,  which  in  turn  was  inversely  proportional 
to  the  meteoroid  flux 


t 

a  p 


oc 


F^V3 


(3) 


The  is  the  meteoroid  flux  having  mass  m^  or  greater.  The  armor  thickness  t 
is  calculated  such  that  particles  smaller  than  m  will  not  penetrate.  Therefore,  F^  is 
also  the  average  puncture  rate  per  unit  of  vulnerable  area.  Use  of  the  general  case 
where 


F^  oc  m 


(4) 


would  have  led  to 


ta  ^  oc  (5) 

Thus  reference  7  has  used  the  special  case  where  jS  =  1.  For  use  where  iS  /  1,  the 
equations  in  reference  7  were  rewritten  to  include  /3.  The  figure  of  merit  r]  is  then  de¬ 
fined  as 


The  relation  between  weight  and  77  then  becomes 

cc  r?(N  a 
N„a„  ® 


(6) 


iV 
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where  is  the  weight  of  the  segmented  prime  surface  radiator.  Reference  7  suggests 
that  for  analysis  of  a  given  vehicle  and  mission  N^a^r  may  be  considered  to  be  con¬ 
stant,  which  makes  77  a  very  useful  figure  of  merit  in  that  it  is  directly  proportional  to 
weight.  Carrying  this  further,  the  ratio  of  segmented  to  unsegmented  weight  can  be  ob¬ 
tained  by 


W* 

s  _  n 


% 


(8) 


where  rj^  is  the  value  obtained  from  equation  (6)  when  the  radiator  is  unsegmented,  so 

that  N  and  N  are  both  unity .  The  values  for  77  presented  in  reference  7  were  for  the 
s 

special  case  where  jS  =  1  and  this  is  the  value  of  the  exponent  used  in  Watson’s  (1956) 
meteoroid  flux  estimate.  Whipple’s  (1961)  meteoroid  flux  estimate  has  a  value  of 
P  =  10/9,  therefore  r?  was  recalculated  using  the  new  value  of  jS;  the  W^/WJ^  was  also 
calculated. 

The  ratio  of  segmented  to  unsegmented  weight  W  /W  ,  which  includes  fins  and  un- 
segmented  manifolds  obtained  by  actually  calculating  the  radiator  weights  for  the  seg¬ 
mented  and  unsegmented  radiator  is  plotted  in  figure  9  (p.  16)  for  the  high  flux  estimate 
and  in  figure  10  (p.  17)  for  the  low  flux  estimate.  Also  plotted  in  figures  9  and  10  is  the 

ratio  W  /W* ,  for  the  prime  surface  radiator,  obtained  from  equation  (8).  Just  as  the 
s  u^ 

survival  probability  P(0)  of  the  unsegmented  radiator  was  set  at  0.  9,  0.  95,  and  then 
0.  99,  the  overall  probability  S  that  out  of  N  segments  as  well  as  the  manifolds 
would  not  be  punctured  was  set  at  0. 9,  0. 95,  and  then  0. 99.  Expressed  mathmatically  for 
N  identical  and  independent  segments. 


N 


S  =  P(0)  X  P(0)j^  =  P(0)  X  y  p''(l 

^s  “  Z-/  n!  (N  -  n): 


n=N. 


P) 


N-n 


(9) 


The  probability  that  the  manifold  would  not  be  punctured  P(0)jj^  and  the  probability  that 

N  of  N  segments  would  remain  unpunctured  P(0)|j  were  equally  weighted.  With  S 
S  iNg 

assigned,  P(0)  and  P(0)j^  are  simply  the  square  root  of  S.  The  individual  segment 
^  s 

probability  of  nonpuncture  p,  which  when  used  in  equation  (9)  yielded  the  desired  value 
of  S,  was  used  when  calculating  tube  armor  thickness. 

Two  generalizations  can  be  made  about  all  of  the  plots  in  figures  9  and  10.  First, 
the  ratio  N  /N  at  which  the  minimum  weight  occurs  is  somewhat  higher  than  the  ex- 
pected  0.75  from  reference  7.  When  using  Whipple's  (1961)  meteoroid  flux  estimate,  the 
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Ratio  of  surviving  to  total  number  of  segments,  N^/N 
(c)  Overall  survival  probability,  0.99. 

Effect  of  number  of  segments  surviving  and  total  number  of  segments  on  segmented  radiator  relative  weight  Whipple's  flux 
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Effect  of  number  of  segments  surviving  and  total  number  of  segments  on  segmented  radiator  relative  weight  Watson's  flux 


minimum  weight  occurred  at  N„/N  of  0.8  or  higher,  and  when  using  Watson's  (1956) 
meteoroid  flux  estimate,  the  minimum  occurred  at  N^/N  closer  to  0.  9.  This  shift  to 
higher  values  is  probably  due  to  the  fact  that  the  manifolds  were  not  segmented.  The 
initial  assumption  that  the  entire  heat  load  be  rejected  from  segments  means  that 
N  -  N  segments  are  redundant  at  the  start  of  the  mission.  Increasing  the  ratio  N  /N 
from  0.  75  to  0.  8  or  0.  9  reduces  the  number  of  redundant  segments,  which  consequently 
reduces  manifold  length  and  weight. 

The  second  generalization  that  can  be  made  is  that  for  any  given  survival  probability, 
number  of  segments,  N^/N  ratio,  and  meteoroid  flux  estimate  the  value  of  de¬ 

creases  as  the  system  power  level  increases;  that  is,  segmenting  is  more  advantageous 
at  higher  power  levels . 

When  figure  9  and  Whipple's  (1961)  flux  estimate  are  used,  it  can  be  seen  that  the 
plot  of  W' /W'  is  a  good  indicator  of  the  relative  weight  for  10  and  20  segments,  but  for 
100  segments  the  level  of  W„/W  is  higher  than  W'/W'.  The  curve  for  W'/W’  ap- 
plies  to  a  radiator  where  100  percent  of  the  weight  is  protective  armor.  When  Whipple's 
(1961)  meteoroid  flux  estimate  is  used  for  the  radiator  with  fins  and  unsegmented  mani¬ 
folds,  the  tube  armor  comprises  only  about  60  percent  of  the  total  radiator  weight.  As 
the  number  of  segments  is  increased,  the  tube  armor  comprises  a  smaller  percentage  of 
the  total  weight,  which  makes  W^/W^  less  applicable  as  a  weight  trend  indicator.  When 
Watson's  (1956)  meteoroid  flux  estimate  is  used  (which  is  considerably  lower  than 
Whipple's,  1961),  the  tube  armor  comprises  a  lesser  percentage  of  the  total  weight  of  the 
radiator  with  fins  and  unsegmented  manifolds.  Thus  W^/W^  as  a  weight  trend  indicator 
is  even  less  applicable  when  Watson's  (1956)  flux  estimate  is  used  for  the  radiator  with 
unsegmented  manifolds  as  can  be  seen  in  figure  10. 

It  is  seen  by  examining  figure  9  that  the  relative  weight  W  /W  for  10,  20,  and 
100  segments  with  Whipple's  (1961)  flux  estimate  is  approximately  0.7,  0.6,  and  0.45, 
respectively,  for  a  design  survival  probability  of  0.  9  and  an  N^/N  ratio  of  0.  8,  is  0.  6, 
0.5,  and  0.4,  respectively,  for  a  design  survival  probability  of  0.95,  and  is  0.45,  0.4, 
and  0.  3,  respectively,  for  a  survival  probability  of  0.  99.  Though  the  weight  of  any  nec¬ 
essary  segmenting  hardware  is  not  included,  substantial  weight  savings  should  be  attain¬ 
able. 

Figure  10  shows  the  results  when  Watson's  (1956)  meteoroid  flux  estimate  was  used. 
For  100  segments,  the  minimum  relative  weight  W  /W  ranged  from  0.  55  to  0.  7  for  a 
design  survival  probability  of  0.  9,  from  0.  48  to  0.  62  for  a  design  survival  probability  of 
0.  95,  and  from  0.  35  to  0.  48  for  a  survival  probability  of  0.  99.  While  the  weight  savings 
of  segmenting  are  not  as  pronounced  as  when  using  the  more  severe  Whipple's  (1961)  flux 
estimate,  segmenting  potentially  can  provide  significant  weight  reduction.  Again  it  should 
be  emphasized  that  the  manifolds  were  not  segmented.  Had  the  manifolds  been  segmented, 
even  lower  relative  weights  would  have  been  obtained. 
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(b)  Overall  survival  probability,  0.99. 

Figure  11.  -  Relative  weight  obtained  when  tube 
length  was  optimized,  compared  to  original 
relative  weight.  Whipple's  flux  (1961). 


A  digression  at  this  point  is  in  order  to  ex¬ 
plain  that  the  segmented  radiator  was  not  opti¬ 
mized  with  respect  to  the  tube  length.  Even 
though  the  tube  length  associated  with  the  mini¬ 
mum  radiator  weight,  for  any  given  size  power 
system,  was  not  the  same  for  each  design  sur¬ 
vival  probability  and  meteoroid  flux,  the  refer¬ 
ence  radiator  was  considered  to  have  tube  lengths 
of  15  feet  for  the  1-megawatt  system,  20  feet  for 
the  2-megawatt  system,  30  feet  for  the 
5-megawatt  system,  and  40  feet  for  the 
10-megawatt  system.  These  lengths  were  se¬ 
lected  from  a  plot  similar  to  figure  5  (p.  9)  but 
for  a  design  survival  probability  of  0.  9.  The 
same  tube  lengths  were  used  for  all  the  segment¬ 
ing  calculations.  The  complete  optimization, 
however,  where  tube  length  was  also  varied  was 
done  for  comparison  by  using  Whipple’s  (1961) 
meteoroid  flux  for  1-  and  10-megawatt-electric 
power  levels  and  for  10  and  100  segments  at  each 
power  level.  This  was  done  for  overall  survival 
probabilities  of  0.  9  and  0.  99.  The  lightest  seg¬ 
mented  weight  at  each  N  /N  ratio  was  used  in 
conjunction  with  the  lightest  unsegmented 
weight.  The  new  values  of  W^/W^  were  then 
superimposed  on  the  original  plots  of  figure  9 
where  S  was  0. 9  (N  =  10,  100)  and  0.  99 
(N  =  10,  100),  and  this  is  shown  in  figure  11. 

The  new  relative  weight  was  the  same  as  for  the 
nonoptimized  plot  for  S  =  0.  9  and  somewhat 
lower  than  the  nonoptimized  plot  for  S  =  0.  99. 
Thus  the  relative  weights  presented  in  figures  9 
and  10  may  be  somewhat  higher  than  if  the  tube 
length  had  been  included  in  the  segmented  radia¬ 
tor  optimization. 

Plotting  the  minimum  weights  for  segmented 
tube  panels  and  unsegmented  manifolds  on  a 
log-log  plot  against  power  level  (fig.  12,  p.  20) 
resulted  in  a  straight  line  as  it  did  for  the  un- 
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(a)  Overall  survival  probability,  0.  9.  (b)  Overall  survival  probability,  0. 95.  (c)  Overall  survival  probability,  0.  99. 

Figure  12.  -  Effect  of  meteoroid  flux  estimate,  number  of  segments,  and  system  power  level  on  segmented  radiator  specific  weight. 

segmented  case,  which  again  indicated  that  the  weight  is  proportional  to  the  power  level 
raised  to  some  exponent.  For  the  unsegmented  case,  the  value  of  the  exponent  b  was 
about  1.4  and  was  influenced  slightly  by  the  survival  probability  and  the  choice  of  meteo¬ 
roid  flux.  When  segmenting  of  the  tube  panels  was  employed,  b  was  reduced  to  lower 
values.  When  Whipple's  (1961)  meteoroid  flux  estimate  was  used,  the  value  of  b  for  all 
three  design  survival  probabilities  was  1.37  for  10  segments  and  1.34  for  100  segments. 
With  Watson's  (1956)  meteoroid  estimate  and  the  same  three  design  survival  probabilities 
the  exponent  b  was  approximately  1.3  for  10  segments  and  1.25  for  100  segments.  The 
trends  are  in  the  expected  direction,  because  increasing  the  number  of  segments  or  de¬ 
creasing  the  severity  of  the  meteoroid  flux  estimate  decreases  the  amount  of  necessary 
protective  armor.  If  no  protective  armor  were  necessary,  the  weight  increase  with 
power  increase  should  be  almost  linear,  which  means  that  b  would  be  unity . 


CONCLUDING  REMARKS 


frhe  salient  results  of  the  study  oriented  principally  toward  radiator  weight  aspects 
of  large  megawatt  size  Rankine  cycle  space  power  systems  are  as  follows  :| 
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I  1.  Radiator  segmenting  can  offer  significant  weight  savings  especially  if  high  survival 
probabilities  are  demanded  and/or  a  high  meteoroid  flux  estimate  is  used.  If  a  radiator 
design  survival  probability  of  0.  99  were  demanded  and  a  severe  meteoroid  flux  such  as 
Whipple's  (1961)  were  used,  100  segments  would  decrease  the  weight  to  less  than  0. 3  of 
the  weight  of  an  unsegmented  radiator.  For  the  same  design  survival  probability  and 
meteoroid  flux  estimate  and  as  few  as  10  segments,  the  weight  of  the  segmented  radiator 
was  still  less  than  half  the  weight  of  the  unsegmented  radiator .  With  a  less  severe  meteo¬ 
roid  flux  estimate  such  as  Watson’s  (1956),  the  weight  savings,  while  not  as  spectacular, 
are  still  large  enough  to  make  segmentation  worthy  of  consideration.  In  this  study  the 
tube  panels  were  segmented  but  not  the  manifolds.  If  the  manifolds  had  been  segmented, 
the  results  would  differ  considerably;  undoubtedly  greater  weight  savings  would  have 
been  possible. 

2.  The  ratio  of  the  design  value  of  surviving  segments  to  total  number  of  segments 
N  /N  that  produced  the  minimum  weight  radiator  was  higher  than  the  expected  value  of 
0. 75.  Values  of  0. 8  to  0.  9  were  foimd  to  yield  minimum  radiator  weight,  but  this  up¬ 
ward  shift  is  probably  due  to  the  fact  that  the  manifolds  were  unsegmented.  Higher  values 
of  N  /N  imply  fewer  redundant  segments,  which  results  in  shorter  and  lighter  manifolds. 
Therefore,  if  the  manifolds  had  been  segmented,  the  optimum  ratio  of  Ng/N  would 


probably  be  closer  to  0.75. 

3.  The  radiator  weight  for  constant  cycle  conditions  was  proportional  to  the  gene¬ 
rated  power  raised  to  some  exponent,  W  cc  p^.  The  exponent  b  was  a  function  of  the 
meteoroid  flux  and  the  desired  survival  probability.  For  the  unsegmented  beryllium 
radiator  at  1400°  F,  the  value  of  b  was  approximately  1.4  when  Whipple’s  (1961)  meteo¬ 
roid  flux  estimate  was  used,  and  just  slightly  less  for  Watson's  (1956)  meteoroid  flux 
estimate.  For  the  segmented  radiator,  the  value  of  b  was  lower  than  for  the  unseg¬ 


mented  radiator 


Lewis  Research  Center, 

National  Aeronautics  and  Space  Administration, 
Cleveland,  Ohio,  October  13,  1965. 
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APPENDIX  A 


manifold  surface  area,  ft 
2 

Ag  surface  area,  ft 

A  area  vulnerable  to  meteoroid 

V  2 

penetration,  ft 

A  ,  internal  surface  area  of  an  indi- 
vidual  tube,  ft^ 

Ug  thin  plate  and  spall  adjustment 

in  eq.  (B7) 

a„  vulnerable  area  in  one  radiator 

^  2 
segment  in  eq.  (7),  ft 

b  ejqjonent  in  eq.  (1) 

C  defined  by  eq.  (H9) 

c  sonic  velocity  in  tube  armor, 

ft/  sec 

Dj  internal  tube  diameter,  ft 

manifold  local  internal  diameter, 
ft 

D  external  tube  diameter,  ft 

o  ’ 

2 

Young's  modulus,  Ibf/in. 

e  element  of  fin 

Ffin  view  factor  from  fin  to  space, 

average 

view  factor  from  tube  to  space, 
average 

F^  meteoroid  flux, 

particles/  (ft^)(yr) 

f  Fanning  friction  factor 

f(Lg)  function  of  generalized  length 

parameter  L 

o 


mass  velocity,  lbm/(ft^)(sec) 

gravitational  constant, 

32.2  ft-lbm/(lbf)(sec^) 

fluid  enthalpy,  Btu/lbm 

Joule's  equivalent, 

778.  2  ft-lbf/Btu 

thermal  conductivity, 
Btu/(hr)(ft)(°F) 

generalized  length  parameter 
defined  by  eq.  (C2) 

jf  distance  from  fin  root  or  distance 

from  manifold  entrance,  ft 

ffin  half  fin  length,  ft 

f  ^  manifold  length,  distance  along 

manifold  axis  between  station 
of  maximum  diameter  and 
station  of  minimum  diameter,  ft 

f  ^  individual  tube  length,  ft 


M 

total  fluid  mass  flow  through  the 
radiator,  lom/sec 

Jr 

local  fluid  mass  flow  in  individual 
manifold,  Ibm/sec 

m 

fluid  mass  flow  through  individual 
tube,  Ibm/sec 

“p 

meteoroid  particle  mass,  g 

N 

number  of  radiator  segments 

number  of  radiator  panels 

% 

Reynolds  number 

SYMBOLS 

G 
Sc 

h 
J 

k 
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Ns 

design  value  for  the  number  of 

temperature  at  midfin, 

surviving  radiator  segments 

To 

fin  root  temperature,  °R 

Nt 

number  of  tubes 

Tw 

tube  internal  wall  temperature,  °R 

n 

constant  in  eq.  (B7) 

Tx 

tube  surface  temperature,  °R 

P 

static  pressure,  psi 

^Tsc 

amount  of  subcooling  of  the  con¬ 

Pe 

system  power  output,  kW  or  MW 

densate  at  exit,  °R 

P(0) 

probability  of  receiving  no 

ta 

armor  thickness,  ft 

meteoroid  puncture  during 

t* 

difference  between  manifold  ex¬ 

P 

mission  duration 

probability  that  individual  seg¬ 

d. 

*fin 

V 

ternal  and  internal  radius,  ft 

fin  thickness,  ft 

ment  will  not  be  punctured  by 

fluid  velocity,  ft/sec 

a  meteoroid  during  mission 

duration 

-r 

3 

material  volume,  eq.  (G3),  ft 

%in 

heat  rejected  per  unit  time  from 

V 

o 

fluid  specific  volume,  ft  /ibm 

individual  fin,  Btu/hr 

V 

assumed  average  meteoroid 

heat  rejected  per  unit  time  from 

velocity,  ft/sec 

manifolds,  Btu/hr 

W 

unsegmented  radiator  weight,  Ibm 

Qtot 

heat  rejected  from  entire  radi¬ 

^a 

tube  armor  weight,  Ibm 

ator,  Btu/hr 

^s 

segmented  radiator  weight,  Ibm 

Qt 

heat  extracted  from  working 
fluid  per  unit  time  per  tube. 

Wu 

unsegmented  reference  radiator 

%in 

Btu/hr 

heat  radiated  from  individual  fin 

X 

weight,  Ibm 

vapor  quality 

per  unit  length  per  unit  time. 

p 

exponent  of  particle  mass  in 

•It 

Btu/  (hr)  (ft) 

heat  radiated  from  tube  external 

y 

eq.  (4) 

coefficient  in  eq.  (B7) 

surface  per  unit  time  per  tube. 

A 

increment 

R 

S 

Btu/hr 

local  radius  of  manifold,  ft 

probability  that  N  segments 
will  not  be  punctured  by  meteo¬ 

e 

V 

\ 

surface  emissivity 

defined  by  eq.  (6) 

value  of  T]  when  N  =  N  =  1 

roids  for  mission  duration 

e 

velocity  exponent  in  eq.  (B7) 

s 

manifold  surface  arc,  eq.  (FI), 

fluid  viscosity,  Ibm/ (ft) (sec) 

ft 
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density,  Ibrn/ft*^ 

Stefan -Boltzman  constant, 

0. 173x10"®  Btu/(hr)(ft^)(°R‘^) 


T  mission  duration  time,  yr 

angles  defined  in  appendix  H 
(p  density  exponent  in  eq.  (B7) 


<P, 


g 


% 


Lockhart -Martinelli  multiplier, 
eq,  (D4) 

two-phase  flow  modulus  for 
turbulent-turbulent  case, 
eq.  (D5) 

ip  angle  in  eq.  (HS) 

n  angle  in  eq.  (H6) 

Subscripts : 
f  liquid 


fg 

g 

It 

M 


refers  to  change  by  evaporation 

vapor 

midfin 

momentum 


m  manifold 

n  number  of  times  an  iterative 

operation  has  been  performed 

p  meteoroid  particle 

T  target  or  protective  armor 

TPF  two-phase  friction 

t  tube 

0  initial 

Superscript: 

'  refers  to  prime  surface  radia¬ 

tor  where  entire  radiator  is 
at  one  constant  temperature 
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APPENDIX  B 


ANALYSIS 

The  calculations  for  this  analysis  were  performed  with  a  computer  program,  but  in¬ 
stead  of  showing  all  the  details  of  the  program,  only  the  basic  equations  that  enable  the 
reader  to  evaluate  the  method  of  analysis  will  be  shown.  Initial  inputs  were  the  radiator 
heat  load,  the  fluid  total  temperature  and  inlet  quality,  the  individual  tube  length,  the 
number  of  panels,  the  desired  tube  pressure  drop,  the  physical  and  thermal  properties  of 
materials,  the  initial  number  of  tubes,  and  a  guess  of  tube  internal  diameter. 

The  description  is  for  the  unsegmented  radiator.  The  differences  between  the  equa¬ 
tions  used  for  the  unsegmented  and  the  segmented  radiators  where  the  tube  banks  were 
segmented  will  be  indicated  at  the  end  of  this  appendix. 

Before  listing  the  equations  a  brief  summary  of  the  approach  used  may  be  of  assist¬ 
ance.  First,  the  individual  tube  mass  flow  or  the  total  flow  through  the  radiator  divided 
by  the  number  of  tubes,  the  individual  tube  length,  and  the  desired  tube  pressure  drop 
were  used  to  find  the  tube  internal  diameter.  Then  knowing  the  total  tube  wall  area  vul¬ 
nerable  to  meteoroid  puncture,  based  on  internal  diameter,  an  armor  thickness  was  cal¬ 
culated,  With  this  first  estimate  for  armor  thickness  the  vulnerable  area  based  on  tube 
outer  diameter  was  used  to  calculate  the  second  estimate  of  armor  thickness.  This  step 
was  repeated  until  the  ratio  of  two  successive  values  was  unity  plus  or  minus  a  1-percent 
tolerance. 

Next  the  heat  radiated  from  the  tube  outer  surface  qj.  was  calculated.  A  surface 
temperature  equal  to  the  internal  wall  temperature,  which  was  assumed  to  be  the  same 
as  the  fluid  total  temperature,  was  used  to  obtain  a  first  estimate  of  the  heat  radiated 
from  the  tubes.  With  the  first  estimate  of  q^.  the  temperature  drop  through  the  armor 
was  calculated.  With  the  new  surface  temperature  q^  was  recalculated.  This  was  re¬ 
peated  until  the  surface  temperature  and  q^  were  compatible. 

Since  q^.  and  the  total  heat  rejection  requirement  were  known,  the  heat  load  that  had 
to  be  rejected  by  additional  finning  was  determined.  Fin  dimensions  were  then  calculated, 
and  with  this  first  estimate  of  fin  length  the  manifold  length  and  surface  area  were  calcu¬ 
lated.  The  armor  thickness  was  recalculated  by  using  the  total  tube  external  surface 
area  plus  manifold  surface  area  as  the  area  vulnerable.  The  manifold  heat  rejection  was 
calculated,  heat  radiated  from  the  tubes  recalculated,  and  the  new  fin  heat  load  used  to 
obtain  new  fin  dimensions.  This  procedure  was  repeated  until  the  fin  geometry  was  such 
that  the  heat  rejected  from  the  tubes,  manifolds,  and  fins  to  space  matched  the  required 
radiator  heat  load.  Then  radiator  weight  was  calculated.  The  total  number  of  tubes  was 
changed  and  the  process  repeated. 

The  equations  used  in  the  analysis  are  now  given.  Mass  flow  through  the  individual 
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tube  is 


m  = 


_M 

N. 


(Bl) 


where  M  is  total  radiator  mass  flow  and  the  number  of  tubes.  Since  tube  mass 
flow  and  tube  length  were  known  and  a  guess  for  tube  internal  diameter  was  made,  the 
tube  pressure  drop  was  calculated  by  using  the  equations  given  in  appendix  D: 


APt  -  -APtpf 


+  P 


M 


(B2) 


which  is  a  sum  of  the  two-phase  friction  and  momentum  pressure  changes.  The  AP^ 
was  tested  to  see  if 

pesired  AP^.  -  0.  05  psi)  <  AP^  <  pesired  AP^.  +  0.  05  psi)  p3) 


K  APj.  did  not  fall  in  the  desired  range,  was  reestimated  by  using 


0.16 


Desired  AP^. 


(B4) 


and  then  AP^.  was  recalculated.  The  exponent  in  (B4)  was  an  arbitrary  choice  based  on 
some  preliminary  calculations. 

When  the  conditions  of  (B3)  were  satisfied,  attention  was  turned  to  the  calculation  of 
the  protective  armor  thickness  required.  For  the  first  estimate  it  was  said  that  the  tube 
outer  diameter  was  equal  to  the  tube  internal  diameter  and  that  manifold  vulnerable  area 
was  zero: 


\,t=VDoNt 

(B5) 

=  Vt+V 

(B6) 

Equation  (B7),  from  reference  11,  for  the  armor  thickness  required  was 
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t 


a 


where 


(B7) 


(The  values  of  these  parameters  were  taken  from  ref.  11  and  used  herein. ) 

The  first  trial  had  been  based  on  =  D^;  therefore,  was  changed  to 

D„  =■  Dj  f  2t^  (B8) 

and  equation  (B7)  was  used  to  get  a  new  t^.  This  was  repeated  until 


^  0.001 

12 


(B9) 


The  heat  radiated  from  each  tube  was 

where  F.  n  =  0-  75  7?  =7?  were  used  as  first  estimates.  Then  the  tube  surface 

tj  0  X  w 

temperature  was  calculated  in  accordance  with  reference  16. 


T 


X 


2irk 


(Bll) 
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When  the  new  value  for  was  used,  was  recalculated  by  using  equation  (BlO). 
This  process  was  repeated  until 


0.999  <  - —  <  1.001 

W„-i 


(B12) 


Then  the  required  fin  heat  load  was  calculated  by  using 


=  1  /^TOT  ^ 


(B13) 


Fin  thicloiess  and  fin  length  were  then  calculated  by  using  equations  (B14)  and  (B15)  with 
a  first  estimate  for 


^fin 


'W 

v^fin 


f(L  )  X  632.  S^CTk  X  2e  X  10^( 

S  VlOOO. 


(B14) 


=  L 

fin  g 


“fm  1 

2crel0® 

^0  f 

Uooo  /  _ 

1/2 


(B15) 


The  development  of  equations  (B14)  and  (B15)  is  shown  in  appendix  C.  With  the  fin  length 
and  known,  the  fin  view  factor  was  calculated  as  shown  in  appendix  H.  Then  fin 
thickness  and  fin  length  were  recalculated  by  using  the  new  view  factor.  This  was  re¬ 
peated  until 


0.999  < 


<  1.001 


Next  the  manifold  surface  area  developed  in  appendix  F  was  calculated: 


(B16) 
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2  ,  t  1  .  2,N  i  /LiL_y  \  ,  (b17) 


For  the  unsegmented  radiator  the  manifold  vulnerable  area  was  included  with  the 
tubing  vulnerable  area  by  returning  to  equation  (B6)  and  recalculating  everything  through 
equation  (B16).  When  the  armor  thickness  based  on  the  combined  tube  and  manifold  area 
had  converged  to  within  ±0. 1  percent  of  the  previous  value  it  was  accepted. 

The  heat  rejected  from  the  manifold  was  calculated  next.  The  view  factor  was  as¬ 
sumed  to  be  0.  85  and  for  a  first  estimate  used  T^^  g  ^  "^w  '^m,f  ^  "^w  "  ^^’^scV^' 


^m,g  ^•^^^^^m,g'^m,g 


(B18) 


Qm,f  "  °-®^'^^^m,f'^m,f 


(B19) 


*^m  ^m,  g  ^m,f 


(B20) 


where  the  subscript  g  refers  to  the  vapor  manifold  and  f  to  the  liquid  manifold.  Then 
the  manifold  external  temperatures  were  calculated  by  treating  the  manifold  as  a  flat 
plate  of  area  g  or  ^  and  thickness  t^; 


rp  _  rp 

m,g=  w--^ 


Q  t 
^m,g  a 


(B21) 


<301.  l‘a 


(B22) 


Then  ^  and  £  were  recalculated  by  using  the  calculated  values  for  T^  ^  and 
Tm  £.  This  was  repeated  until 


0. 999  < - ^  <  1.  001 


(B23) 


The  heat  to  be  rejected  by  the  fin  was  recalculated; 
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(B24) 


Qtot  ■ 


ItNt 


The  F|.  was  calculated  as  shown  in  appendix  H.  With  the  new  value  for  the  pro¬ 
gram  returned  to  equation  (B14)  to  calculate  a  new  Everything  from  equation  (B14) 

to  (B24)  was  repeated  until 


0.999  < 


<  1.001 


(B25) 


Then  the  total  vulnerable  area  was  calculated  and  checked  against  the  previous  value  to 
see  that 


0.999  < 


(B26) 


If  the  condition  of  equation  (B26)  was  not  met,  the  program  returned  to  equation  (B6). 
When  the  ratio  fell  within  the  range  specified  in  equation  (B26),  the  manifold  pressure 
drops  and  radiator  weight  were  calculated. 

The  pressure  drops  in  the  manifolds  were  calculated  by  the  equations  developed  in 
appendix  E.  The  fluid  velocities  and  were  assumed.  The  actual  values  used 
are  listed  in  the  assumptions  in  the  text.  The  working  fluid  properties  were  obtained 
from  reference  15: 


AP 


^g,  0  g  g  ni 


m,g 


144  g, 


2M 


1/2 


0.  0028  + 


0.3676 


2Mp„V, 


■nil, 


0. 16 


AP 


m 


m,f 


144  g, 


4M 


.1/2 


,N.^irPfV, 


0.  0028  + 


0. 3676 


'4MpjVj 

N 


0. 16 


(B27) 


(B28) 
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(B29) 

(B30) 

(B31) 

(B32) 

(B33) 

(B34) 


Modifications  for  Segmenting 

Segmenting  was  employed  to  reduce  the  amount  of  protective  armor  necessary  to 
maintain  the  same  radiator  survival  probability  as  for  the  unsegmented  radiator.  Re¬ 
dundant  radiator  sections  were  carried  so  that  if  an  active  radiator  section  were  punctured 
by  a  meteoroid  the  damaged  section  would  be  isolated  and  one  of  the  redundant  sections 
activated.  Thus  the  radiator  with  N  segments  was  designed  to  reject  the  entire  heat 

load  from  N  segments,  N  being  the  design  value  for  the  number  of  segments  which 
s  s 

would  remain  unpunctured  at  the  end  of  the  mission.  This  requires  minor  modifications 
of  the  equations  used  for  the  unsegmented  radiator.  First  of  all,  since  only  seg¬ 
ments  out  of  a  total  of  N  segments  are  active,  the  number  of  active  tubes  is  the  frac¬ 
tion  N  /N  times  the  total  number  of  tubes.  Thus  in  equations  (Bl),  (B13),  and  (B24) 
s 

where  appears  it  is  replaced  by  N^(Ng/N). 

The  other,  and  most  important  difference,  is  the  procedure  for  calculating  the  armor 
thicknesses.  The  armor  thicloiess  for  the  manifold  was  calculated  separately  from  the 
tube  armor  thickness.  This  was  necessary  because  the  manifolds  were  not  segmented. 


only  the  tube  banks.  The  overall  survival  probability  S  that  the  manifolds  and  out 
of  N  independent  segments  would  remain  unpunctured  is 


N 


s  =  P{0)jnP(0)N  " 


m 


Z 


n=N 


N, 


(N  -  n);n; 


p"(i  -  p)"-" 


(B35) 


The  value  of  S  was  assigned  and  P(0)^,  the  probability  that  the  manifolds  would  not  be 

punctured,  and  P(0)j^  ,  the  probability  that  of  N  segments  would  remain  unpunc- 

s 

tured,  were  weighted  equally,  or 


P(0)j„  =  P(0)n  (B36) 

S 

A  value  of  p,  the  probability  that  any  individual  segment  would  not  be  punctured,  was 
found  so  that 


mu 


s 


Z 


n; 

(N  -  n).'n.' 


p^l  -  p)N-''=V5 


(B37) 


The  tube  vulnerable  area  of  an  individual  segment  is 


a 


V 


^t 

N  ‘  ° 


(B38) 


Then  the  tube  armor  thickness  would  be  calculated  by  using  equation  (B7)  with  this  modi¬ 
fication.  The  term  [^o'A^Ty'-ln  was  replaced  by  the  term  ^aa^ry-ln  pj 

Similarly  the  manifold  armor  thickness  was  calculated  by  using  equation  (B7)  with  the 
term  P(0)J^^^^  replaced  by  j^ofAj^T/'ln  Equation  (B6)  is  de¬ 

leted  when  calculating  the  segmented  radiator.  Essentially  these  were  the  only  impor¬ 
tant  differences  between  the  segmented  and  unsegmented  calculation  procedures. 
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APPENDIX  C 


FIN  LENGTH  AND  THICKNESS 

The  equations  for  calculating  half  fin  length  and  fin  thickness,  with  fin 

width  fj.  held  constant  were  taken  from  reference  13.  One  modification,  adding  a  fin 

view  factor,  was  made  to  account  for  the  fact  that  the  fin 
did  not  have  an  unrestricted  view  of  space,  see  sketch  (a). 


Substituting  equation  (C2)  into  (Cl)  and  simplifying  yield 


Since  L  has  been  chosen  and  the  equivalent  sink  temperature  was  assumed  to  be  zero, 
T£  can  be  obtained  from  figure  8  of  reference  13.  With  known  and  l(Lg)  obtained 
from  equation  (39)  of  reference  13,  *fin  can  be  found  from  equation  (C3).  It  remains 
then  to  solve  for  in  equation  (C2). 
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APPENDIX  D 


TWO-PHASE  PRESSURE  DROP  IN  A  CONSTANT  DIAMETER  TUBE 

The  flow  regime  assumed  was  that  of  a  homogeneous  mixture  with  both  the  liquid  and 
vapor  traveling  at  the  same  velocity.  To  obtain  two-phase  pressure  drop,  the  single¬ 
phase  multiplier  of  reference  17  was  used. 

The  tube  was  divided  into  increments  and  the  two-phase  friction  static  pressure  drop 
for  each  increment  was  calculated  by  using  the  local  vapor  velocity.  Because  the  con¬ 
densing  process  is  essentially  isothermal  the  wall  heat  flux  was  deemed  to  be  constant. 

The  energy  equation  as  applied  to  the  increment  is 


Qt 

4—  A£ 


V 


^n-1  + 


n-1 


w 


2gcJ 


3600  GD. 


(Dl) 


For  the  initial  calculation  the  velocity  term  was  neglected  but  then  considered  after  the 
local  vapor  quality 


n 


was  found.  The  average  quality  for  the  increment  is  then 


Vn-l 


^n-l+^n 


(D2) 


(D3) 


The  use  of  the  average  quality  will  be  evident  after  showing  the  pressure  drop  equation. 
From  reference  17 


TPF 


where  is  a  function  of 


X- 


For  the  assumed  turbulent -turbulent  case, 


(D4) 
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(D5) 


The  fluid  density  p  and  viscosity  p  were  evaluated  at  inlet  total  temperature  and  pres¬ 
sure,  The  weight  flow  of  the  liquid  and  the  weight  flow  of  the  vapor  ihg  can  be  ex¬ 
pressed  as  a  function  of  the  vapor  quality 

riij  =  (1  -  x)m  (D6) 


and 


m 


g 


=  xm 


Thus, 


Wn,  n 


1 


A  curve  fit  for  the  plot  of  (p^  against  Xtt  obtain  cp^.  Then 


2P, 


At 


144g^D. 


p7) 


(D8) 


(D9) 


where 


n,n-l 


=  G^,n-l''g 


and 


fn.n-l  =  0.0014  + 


0. 125 


"i^g^)n.n-l 


0.32 


(DIO) 


(Dll) 
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Equation  (Dll)  is  the  empirical  relation  of  Koo  (ref.  16)  and  is  valid  for  the  range  of  Np^ 
from  3  000  to  3  000  000.  This  was  chosen  instead  of  the  Blasius  friction  factor,  recom- 
mended  in  reference  17,  because  it  was  applicable  over  a  wider  range  of  Reynolds  num¬ 
bers.  Then  h^^  was  corrected  for  the  change  in  velocity  across  the  increment  and  the 
pressure  drop  for  the  next  increment  found.  The  friction  two-phase  pressure  drop  was 
added  to  the  momentum  pressure  rise: 


AP. 


-APtpp  +  APjyf  ^ 


2('Pgf 


lAi>t 


144  g^D. 


144g, 


(D12) 


n=l 
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APPENDIX  E 


MANIFOLD  PRESSURE  DROP 

It  was  approximated  that  the  fluid  mass  flow  in  the  manifold  decreased  linearly  with 
manifold  length.  The  mass  flow  at  any  station  H  is  then 


where  jf  denotes  the  local  fluid  mass  flow  in  the  individual  manifold.  For  the  vapor 
manifold, 


and  for  the  liquid  manifold 


^0- 


=  2R 


Secondly,  it  was  stipulated  that  the  manifold  diameter  would  be  varied  to  maintain  a  con¬ 
stant  fluid  velocity  for  the  decreasing  mass  flow. 

The  two-phase  friction  pressure  drop  in  the  vapor  manifold  is 


2(plipvl 

g  g 


df/Tpp  144g^^j 


When  Koo's  empirical  equation  (ref.  16),  is  used 


f  =  0. 0014  + 


0. 125 


N' 
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which  is  valid  for  3  000  <  Nj^  <  3  000  000  and 


'R 


g 


(E4) 


Substituting  equations  (E3)  and  (E4)  into  (E2)  yields 


0  .2^  ir2 

144g^®i 


0. 0014  + 


0. 125 


^g^g 


0.  32 


M, 


^9' 


g 


Mass  flow  is 

4  i^g  g 


{E5) 


(E6) 


Rearranging  (E6)  and  substituting  (El)  yield 


4^ 


1/2 


(E7) 


Substituting  equation  (E7)  into  (E5)  and  assuming  that  the  change  in  vapor  quality  is  small 
enabling  p  p  and  (p  in  addition  to  V  to  be  treated  as  constants  yield  equation  (E 5) 

o  o  o  o 

as 
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Integration  of  equation  (E8)  for  the  vapor  manifold  friction  pressure  drop  and  substitution 
for  J(q  yield 


AP 


ni,g 


144g, 


2M 


\i/2 


0.  0028  + 


0. 3676 


For  the  liquid  manifold,  the  friction  drop  is 


2MpgVg 


■nil, 


,0. 16 


AP 


2pjVp 


m 


m.f 


144g, 


4M 


,1/2 


^^rNpPfVf 


0.  0028  + 


0. 3676 


4MpjVj 

V"? 


0.  16 


(E9) 


(ElO) 
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APPENDIX  F 


MANIFOLD  SURFACE  AREA 


The  manifold  external  surface  area  must  be  calculated  and  included  as  part  of  the 
area  vulnerable  to  meteoroid  puncture.  Also,  it  is  required  to  calculate  the  quantity  of 
heat  radiated  to  space  directly  from  the  manifold. 

Using  the  notation  shown  in  sketch  (b)  the  external  surface  area  is 


In  equation  (E7)  it  was  shown  that 


(FI) 


(F2) 


(F3) 
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and 


(l  -  — 

4irpVl^  V  ^ra/ 


Substituting  equations  (F2)  and  (F4)  yields  equation  (FI)  as 


Integration  of  equation  (F5)  yields 


This  cumbersome  expression  was  simplified  after  closer  examination.^  First,  >  1, 
therefore  »  1  and  it  was  found  that  J( ^jfynV  «  1;  therefore,  • 

The  troublesome  term  in  equation  (F6)  was  the  inverse  hyperbolic  tangent.  As 


•^0 
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the 


tanh 


OO 


4?rpV^tn 


However,  looking  at  the  expression  practically,  if 


AirpVit 


>  2x10' 


m 


then 


tanh’ 


K  5.7 


^0 


^irpYi. 


To  simplify  then,  the  insignificant  terms  in  equation  (F6)  were  deleted  yielding 


Ag  ^  2'!t< 


(«m) 


3/2 


12£ 


m 


(F7) 


or 


Ag  ^  27r£^ 


i  (rfo  r  . 

3  VpV/  ^ 


(F8) 


The  manifold  surface  areas  then  are 


M 


1/2 


3  V2.PgVg 


+  t. 


(F9) 
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Am,f  = 


2  /  M  V 


1/2 


3  N 


p^PfVj 


+  t. 


(FIO) 


'^ni,g  *  ^m,f 


(Fll) 
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APPENDIX  G 


MANIFOLD  WEIGHT 

As  in  appendix  E  this  is  the  development  of  the  general  case  for  the  individual  mani¬ 
fold  using  Ji  to  denote  local  fluid  mass  flow  in  the  individual  manifold.  The  sketch  in 
appendix  F  is  also  applicable  here. 

The  volume  of  material  is 


"V  = 


de 


(Gl) 


assuming  t* 


t  .  Simplifying  and  substituting  equation  (E7)  into  (Gl)  yield 
a 


=  JT 


2t. 


since  the  weight  is 


W 


m 


Prj.'f' 


(G2) 


(G3) 


Substituting  equation  (G2)  into  (G3)  and  integrating  yield 


m 


/ ^  ^ 

sXpnV/ 


When  is  replaced  with  the  proper  equivalent,  the  manifold  weights  are 


(G4) 


Wm,g  =  2PT’^y 


m 


4  /  M 


1/2 


3 


+  t. 


and 


^m,f  "  ^p^T^a^ 


m 


4  /  M 


:/r  \l/2 


3  ^NpPfTTVf 


+ 1. 


(G5) 


(G6) 
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APPENDIX  H 


ANGLE  FACTORS 


The  angle  factor  for  the  heat  emission  from  a  differential  element  of  the  fin  to  the 
tube  has  been  derived  as  a  general  case  in  reference  18.  Given  two  surfaces,  one  of 
differential  size  and  the  other  a  cylindrical  surface  of  arbitrary  profile,  whose  generating 
lines  are  parallel  to  the  first  surface,  shown  in  sketch  (c),  it  was  shown  that  the  angle 
factor  of  the  cylindrical  surface  A2  from  the  differential  element  dA^^  is 


^dAi,2  = 


i-  (sin  -  sin  $”) 
2 


(HI) 


For  the  fin  and  tube  geometry  of  the  radiator  =  7r/2.  The  angle  factors  from  the  dif¬ 
ferential  element  of  the  fin  de  to  the  tube  surfaces  designated  1  and  2  in  sketch  (d)  are 


^de,  1  =  ^  (1  -  cos  $1)  (H2) 

and 

^de,2  ^2) 

For  this  study,  the  interradiation  between  the  fin 
and  the  tube  was  neglected  and  only  the  radiation  to 
space  considered.  Thus, 

F  ,  =  1  -  F  1  -  F  o  =  -  (00s  #1  +  cos  #,) 

de, space  de,  1  de,  2  g  1  ^ 


(H4) 
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The  angle  factors  from  f  =  0  to  £  =  were  found  in  increments  and  a  mean  value 
was  calculated  by  using  the  trapezoidal  rule: 


n=£fin/Ai 


de,  space)j^_ j  "*■  de,  space)^^ 


I 


fin 


(H5) 


The  angle  factor  from  any  point  on  the  tube  to  the  fin  and  adjacent  tube  is  found  similarly. 
Here  the  differential  element  (see  sketch  (e))  lies  on  the  tube  and  the  cylindrical  surface 
is  the  fin  surface  plus  part  of  the  adjacent  tube: 

1  =  J  (1  -  cos  n)  (H6) 


For  emission  directly  to  space, 


-  ^  ■  ^dib.  1  -  7  (1  +  cos  S2) 


d:/),  space  d:/),  ^  2 


(H7) 


The  relation  between  and  tp  is  shown  in  reference  19  to  be 


cos  Q 


_  C  sin  :/)  -  2C  cos  ip  -  {C  cos  ip  -  1) 


(H8) 


-  2C  cos  :/)  +  1 


■ 

■ 

B 

m 

b 

i 

a 

1 

(e) 
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The  mean  value  of  the  tube  angle  factor  for  0<il><  ir/2  was  obtained  numerically  by 
using  the  trapezoidal  rule: 
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aeronautical  and  space  activities  of  the  United  States  shall  be 
conducted  so  as  to  contribute  ...  to  the  expansion  of  human  knowl- 
edge  of  phenomena  in  the  atmosphere  and  space.  The  Administration 
shall  provide  for  the  widest  practicable  and  appropriate  dissemination 
of  information  concerning  its  activities  and  the  results  thereof.^* 

— ^National  Aeronautics  and  Space  Act  of  1958 
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